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Abstract. Both observational and theoretical rates of binary neutron star coales- 
cence give low prospects for detection of a single event by the initial LIGO/VIRGO 
interferometers. However, by utilizing at the best all the a priori information on the 
expected signal, a positive detection can be achieved. This relies on the hypothesis 
that 7-ray bursts are the electromagnetic signature of neutron star coalescences. The 
information about the direction of the source can then be used to add in phase the 
signals from different detectors in order (i) to increase the signal-to-noise ratio and 
(ii) to make the noise more Gaussian. Besides, the information about the time of 
arrival can be used to drastically decrease the observation time and thereby the false 
alarm rate. Moreover the fluence of the 7-ray emission gives some information about 
the amplitude of the gravitational signal. One can then add the signals from ~ 10 4 
observation boxes (~ number of 7-ray bursts during 10 years) to yield a positive de- 
tection. Such a detection, based on the Maximum a Posteriori Probability Criterium, 
is a minimal one, in the sense that no information on the position and time of the 
events, nor on any parameter of the model, is collected. The advantage is that this 
detection requires an improvement of the detector sensitivity by a factor of only ~ 1.5 
with respect to the initial LIGO/VIRGO interferometers, and that, if positive, it will 
confirm the 7-ray burst model. 

A widely spread out gamma-ray burst model is related with the coalescence 
of two neutron stars (N.S.) at cosmological distance ||. In this model, the 
gravitational energy liberated during the disruption of the less massive star of a 
binary system of N.S. is transformed into electromagnetic energy and radiated 
in the X,7 range jjj. 

No direct evidence of the validity of this model exits until now. Only energy 
budget considerations and an estimated rate of coalescences of N.S. computed 
in the frame of stellar evolution theory, confort the idea that at least a large 
fraction of 7-ray bursts are generated by the above mechanism Bj . 

In this communication, we want to show that the class of a typical noise 
of 10~ 23 /\/Hz G.W. detectors is sensitive enough to detect, after a few years 
of observation, the associated G.W. emission of the coalescence. The strategy 
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Figure 1: Fraction of the sky covered when phasing the signal of different detec- 
tors (LI = LIGO 1, L2 = LIGO 2, V = VIRGO), in the case where the direction 
of the source is known with an accuracy of 5°. A portion of the sky is said to be 
covered if the phase lag between different detectors resulting from the incurate 
knowledge of the source position is lower than 7r/2. 

proposed here consists in using at the best all the a priori information and 
hypothesis of the model and to "summ" the signals of many events. 

Historically, coalescing N.S. were considered as the most promising source 
of detectable G.W. A detection rate of few events/year was estimated in an (to 
much) optimistic case (for a review on the subject see [Q). Since then, there 
is widely spread out consensus that the coalescing rate of N.S. is one per year 
in a sphere of 200 Mpc of radius. 7-ray burst have a repetition rate 10 times 
lower, and a focusing mechanism must be invocated in order to explain this 
discrepancy [||. 

The predicted sensitivity of the first VIRGO detector will allow us to detect 
a coalescence of N.S. at a distance of 27 Mpc with a signal/noise (S/N) ratio 
of 7 (the figure 7 was choosen in view to have one false alarm per year 
under the hypothesis that the noise is Gaussian). Therefore the detection rate 
of coalescing N.S. will result to be 0.5 event per century (!) or 10 times lower if 
no beaming mechanism exist in the 7-ray burst emission. 

In view to increase the detection rate it is convenient to add in phase the 
signal of different detectors. In fact if N is the number of detectors, the S/N 
increases as the y/N and the detection rate as N 3 / 2 . Moreover, adding different 
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Figure 2: Same as Fig. 1 when the direction of the source is known with an 
accuracy of 1° instead of 5°. 

signals has the big advantage to make the signal more Gaussian (Central limit 
theorem !) and therefore to decrease the false alarm rate. 

Adding in phase the outputs of different detectors is easy if the position of 
the source in the sky is known with a sufficient accuracy. This is the case of an 
optical or radio detection of the source (pulsars or supernovae). 

For 7-ray bursts, the error box of the recent BATSE experiment aboard the 
Compton-GRO satellite is about 5° x 5°. This accuracy in the determination of 
the position of the source in the sky will allow the phasing of the signals from 
3 detectors (VIRGO, LIGO 1 and LIGO 2) for signals whose frequencies are 
lower than 200 Hz (See Fig 1). Figure 2 shows that reducing the error boxes 
by a factor of 5 will be enough to phase the output of 3 detectors. It is quite 
possible that the new generation of X and 7-ray detectors will improve the 
localization accuracy. In our opinion, the Gravitational Wave Community can 
play an important role in convincing the 7-ray Community of the importance 
of improving the accuracy of positioning the sources. 

The hypothesis that 7-ray bursts are generated by coalescing N.S. allows 
us to suppose that the frequency of the gravitational radiation sweeps the fre- 
quency window of the detectors a few seconds before the begining of the 7-ray 
emission; moreover, we suppose that the strength of the gravitational signal 
and the fluence of the electromagnetic one are correlated with the distance of 
the source. These quite natural hypotheses will be considered as a priori in- 
formations and will be used in computing the a posteriori probability density 



of detection of the signal. Other a priori informations will be enumerated (and 
used) in what follows. 

The Maximum a Posteriori Probability Criterium (M.P.P.C.) will be the 
tool used in view to utilize at the best all the a prior information given by the 
model. Let us consider a simple example to show how the M.P.P.C. works. 

Consider a sinusoidal signal of length T: y(t) — ax cos(ujt + cf>) in presence 
of a white Gaussian noise n{t). The sampled signal is: 

Ui = a x cos{uoti + (j>) + rii, i=l,2...,N, (1) 

where N is the number of sampling points: N — T x v s , v s being the sampling 
frequency. The density probability of having N yi values is 

1 - 

P{yi,y2,...y N \a,4>) = - 7=^ exp (- V(y; - acos(ujt t + cj>)) 2 /2<j 2 ) . (2) 

(vV 27T ) i=i 



Let us suppose that the phase 4> is known. In this case the maximum likelihood 
criterium allows us to estimate the amplitude of the signal a. 



0, from which we obtain the well known result a = 2/N^2,f =1 yiCOs{u>t 



The best estimation of a is given by solving the equation d a P(yi, 2/2, ■■■Un\o., 
'rom which we obtain the wel 
and the Signal/Noise ratio reads 

S/N = aT/B , (3) 

where T is the time length of the signal (sec.) and B is the instrumental noise 
per \/Hz- 

Consider now the little more complicated case, in which the phase <f> is not 
known. We have two possible strategies: the first one consists in estimat- 
ing simultaneously a and cj) by solving the system d a P(yi, 2/2, ...yjv|a, </>) = 0, 
d ( pP(yi, 2/2, ■■■VN\a, 4>) — 0. In this way we obtain one value for a and a value 
for <f). This gives too much information because we do not need to know </>. The 
alternative strategy consists in maximizing the posteriori probability density de- 
fined by P{<t>)P{yi , 2/2, ■■■UN\a, 4>)d<j) where P((j>) is the a priori probability 
density of the phase. Because we have no information on the phase (except 
that it spans the interval 0, 2tt[) we apply the so-called Equipartition of the 
Ignorance Principle J?]] , jj] , [ | to set P((f>) — const. — l/2ir. After integration, 
the amplitude a is found by the maximum of the function exp a x Iq (a) where 
I is the Bessel function of order 0] . The asymptotic value of a for S/N — > 00 
is a = y/c 2 + s 2 where c = l/NJ2 i=1 yi cos(wii) and s = V% sin(u;fj), 

i.e. the value of a is given by the power spectrum of the signal. The S/N is 



S/N = aT/(V2 x B) , 



(4) 



i.e. l/\/2 times worse than in the previous case in which the phase was known. 
The factor y/2 is the price to pay for our ignorance. If more then one "box" 
exist, say J boxes, then the M.P.P.C. tell us that we have to maximize the a 
posteriori probability given by the product of the density probability of each 
box: 

Pi(yi,V2, ■•4/Jv|a,</0- p 2(s/2,3/2) -VN\a,<f>)-Pj(yi,y2, -VN\a,<f>) (5) 

It is easy to show that, roughly speaking, the signal-to-noise ratio is increased 
by a factor J 1 / 4 . 

In the real case, i.e. the detection of the signal from coalescing N.S., the 
probability density distribution P{yi,y2, ••■2/iv|a, <t>, J7ii, ma, i.-) depends on more 
variables than a and <j), namely the masses mi and mi of the two N.S., the 
angle i between the orbital plane and the line of sight, etc... The a priori 
informations that we have about these parameters must be used (for example, 
we know that the masses of the two N.S. cannot be arbitrary, but lie between, 
say, Mi = 0.2M© and M 2 = 3M Q , and so on). Therefore, the a posteriori 
probability density reads 



27r pM 2 pM 2 i-ir/2 
J Mi J Mi JO 



P(m 1 )P(m 2 )P(i)P(4>) ... x 
xP(yi,y 2 , ...j/jv|o, 4>, mi,m 2 , i..) d<j)dmi dm 2 dfy 



Under the conservative hypothesis that gravitational radiation observations 
must start, say, one minute before the electromagnetic counterpart of the 7-ray 
burst (the time at which the two N.S. merge is not precisely known either for 
observational reasons and for the lack of a reliable model) , the total observation 
time is decreased by a factor of ~ 2 x 10~ 3 (1 minute/ (number of minutes per 
day) x (number of 7-ray bursts per day)). Therefore the false alarm probability 
is reduced of the same factor (vice-versa S/N is reduced to 6 for the same false 
alarm probability). Note that the idea of using the information given by the 
timing and direction of the 7-ray bursts to increase the detection rate has been 
already developed by Kochanek & Piran ||. 

With five detectors (2 + 2 for LIGO 1 and LIGO 2 and 1 for VIRGO) the 
distance at which a N.S. coalescence can be detected with a S/N of 7 is increased 
by factor y/5, i.e. 60 Mpc or 420 Mpc for S/N = 1. If the coalescence happens 
at a distance of, say, 4.2 Gpc, the S/N ratio will be 0.1. Taking into account 
that in 10 years of observation time, 10 4 7-ray bursts will be observed (The 
7-ray bursts rate is about 3 per day), we shall dispose of 10 4 "boxes" with a 
gravitational S/N ratio of ~ 0.1 in each box. 

By "adding" the signal of the 10 4 boxes, as explained in the above example, 
the S/N ratio will be increased by factor of 10 (= (10 4 ) 1 / 4 ), i.e. S/N = 1. Note 
that this is a pessimistic estimation: in fact not all 7-ray bursts are situated at 



4.2 Gpc. Under the hypothesis of an uniform spatial distribution, the average 
7-burst distance is reduced by a factor v3 and therefore the signal-to-noise ratio 
increases by a factor of-\/3: S/N = 1.73. This is not the end of the story: in fact 
we have not used yet the hypothesis that the intensity of the gravitational signal 
is correlated with the electromagnetic one (more precisely with its fluence) . 

The operation of "adding" the different boxes must be done by weighing 
the boxes with a weight proportional to the fluence of the 7-ray bursts. The 
final result is S/N — 5.5. This is not enough to have an acceptable false alarm 
probability; in fact a S/N > 7 is required to have a false alarm probability 
< 1(T 3 . 

It appears that an improvement by a factor > 1.3 of the sensibility of the 
detectors will be enough to test the hypothesis that 7-ray bursts are generated 
by coalescing N.S. 

The above results must be considered as preliminary ones. In fact we have 
supposed that the detectors are aligned (optimistic hypothesis); in real calcu- 
lation the directionality of the detectors should be taken into account and used 
to weigh the signal of each box. The S/N ratio should be improved a little bit. 

CONCLUSIONS 

There is a general agreement that N.S. coalescence rate is ~ 1 per year 
within a sphere of 200 Mpc. Taking into account the incertitude of the above 
estimation (a factor of 2 - 3, L. Yungelson private communication), this value is 
in a good agreement with the observed rate of 7-ray bursts in the same volume 
(0.1 per year). A possible beaming of a factor 10 is plausible and is invocated 
in view to eliminate the discrepancy between the theoretical model and the 
observational data. 

A more important beaming of the 7-ray emission during the merging of two 
N.S. would increase the coalescence rate. It is however quite difficult to imagine 
which physical mechanism would decrease the isotropy of the electromagnetic 
radiation: the velocity of the two merging N.S. is indeed about half the velocity 
of light and therefore the aberration effects are not strong. 

Recently, Wilson et al. jl0) have found by numerical simulations that the 
more massive of the two N.S. can collapse and form a black hole before the 
merging if its mass is close to the critical one. Even in this unlikely particular 
case the electromagnetic radiation is not suppressed: in fact the less massive 
star will be destroyed by the tidal forces of the just born black hole. 

An improvement of the detectors sensitivity by a factor of at list 10 is re- 
quired in order to detect a few individual events per year with a reasonable 
false alarm probability. Nobody knows when and how the thermal noise of the 
mirrors will be reduced by such an important factor (one of us (S.B.) is indebted 
to Prof. A. Giazotto for helpful discussions on this point). These (quite pes- 



simistic) conclusions hold only for the coalescence of N.S. Let us recall that a 
coalescence of massive black holes (100 M ) can be detected up to <~ 3 Gpc. 

In this paper we have showed that by phasing the signals of different detec- 
tors and by using all the a priori informations given by the 7-ray burst model, 
a positive detection can be achieved in a few years of observation time with a 
little improvement (less than a factor of 2) of the first generation of the VIRGO 
class detectors. The price to pay is the loss of information: at the end of the 
observation only a positive detection of gravitational radiation will be achieved 
and a model for 7-ray burst confirmed, but the time and the position of the 
events will be lost. 
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